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An analysis of the analytical solution of the Schrodinger equation (which is a second order dif- 
ferential equation) for H2 shows that the second linear independent solution of this equation is a 
square integrable function and therefore the ground state total wave function is a linear combination 
of two linear independent wave functions of different space symmetry: cylindrical and spherical. The 
wave function of cylindrical symmetry is well known. It has maxima at the positions of nuclei. The 
wave function of spherical symmetry and the corresponding spherical electron distribution, which 
exists at _R 7^ and locates at the middle of the bond, represents a quasiatom of electron density 
of non-nuclear united atom. In the light of the new result the qualitative behavior of the ground 
state wave function and the electron density of H2 has been reinvestigated. It is shown analytically 
that a transformation of the total molecular wave function with two maxima to that one with one 
maximum passes through a flat wave function. The presented three-dimension figures of the electron 
density visualize the spherical component of the total wave function and its transformation with 
increasing internuclear separation. 

PACS numbers: 31.10.-l-z, 03.65.-w 
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The property and behavior of the ground state wave 
function of H2 molecular ion is very important for quan- 
tum mechanics because it establishes our general repre- 
sentation on the quantum interaction in molecules. For 
this reason the ground state wave function and electron 
density of H2 ion was investigated in many publications 
by using different approaches. However, an analysis of 
these results shows that there is a direct contradiction 
between them. This lead us to non-unique qualitative 
representation on a behavior of the ground state wave 
function. 

Thus, from one hand side, from asymptotic proper- 
ties of the total wave function follow that at the limit 
of a united atom the ground state wave function of H2 
transforms to the wave function of He^ . Therefore, af- 
ter some small R the H2 wave function must have only 
one maximum at the middle of the bond. This is a basis 
of perturbation theory developed by Bethe for a calcula- 
tion of the H2 total energy at short R 1] . In this the- 
ory the H2 wave function is approximated by the He^ 
wave function located at the middle of the bond. The 
investigation of Bethe's perturbation theory shows that 
it works well up to i? = 0.05 a.u. 2]. Other example, 
which displays an importance of a correct asymptotic of 
the total molecular wave function in the limit of a united 
atom was given in Ref. 3]. In this work it was shown 
that a representation of the molecular wave function as a 
linear combination of united and separated atoms wave 
functions results in significant improvement of the H2 
variational total energy calculated at short R. 

From other hand side, from a solution of the 
Schrodinger equation in Born-Oppenheimer approxima- 



tion in alliptic coordinates for if^ given by Burrau [J 
follows that the solution has two maxima at the position 
of two nuclei. 

To clarify this question and to show that the second 
linear independent solution of the Schrodinger equation 
is square integrable function let us consider first the H2 
problem qualitatively. In the Born-Oppenheimer approx- 
imation the Hamiltonian of iJ^ can be written as 
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where r is the vector of the electron; R —\ Ra — Rb \ is 
the internuclear distance; Ra and Rb are the vectors of 
nuclei A and B; and the electron and nuclei charges are 
equal to one. A parametric dependence of H{R) on R 
shows that at the limit i? -^ Hamiltonian (QJ reduces 
to the Hamiltonian of the iJe+ 
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while the H2 total wave function transforms to the total 
wave function of the He^ ion. At the limit i? — > 00 there 
is an equally probability that the electron can be located 
near nucleus A 01 B and thus forms H atom at the point 
A or B. Therefore, Hamiltonian |^ transforms to the 
Hamiltonian of the H atom 

fi2 1 

H{R)^-^V^-^^— 

2m I r — Ra \ 

The total molecular wave function in this case reduces to 
the wave function of the H atom. 

Summarize the asymptotics of Hamiltonian (^ and fol- 
low to the superposition principle we can conclude that 



the total wave function of 7J^ molecular ion can be pre- 
sented as a linear combination of three atomic wave func- 
tions 

*(i?)=CA(il')*A + Cs(i?)*B+Ct/(il')^C/ (2) 

where ^^i and ^>b are the wave functions of H atoms 
and ^[/ is the total wave function of He~^ ion. In this 
expression each wave function '5^, "^b, and ^[/ is located 
on a separated center. Therefore, we can see that con- 
sistent quantum mechanical consideration of ifj prob- 
lem results in a three-center problem with corresponding 
three-center representation of the total molecular wave 
function. 

The use of representation (O for the total molecular 
wave function in calculations of H2 total energy with 
high precision at short R was advocated in Ref. j3|]. At 
equilibrium distances the use two-center expansion in- 
stead three-center expansion Q results in a small gap be- 
tween the Hartree-Fock total energies of H2 calculated by 
numerical methods and in conventional two-center LCAO 
approximation. In Ref. [3 was shown that the use of 
expansion |(2Jl permits to remove this gap. These two ex- 
amples together with the Bethe's perturbation theory yj 
show that the expansion (O is correct. 

One of a particular interest of the present study con- 
sists in study of a transformation of the molecular wave 
function and electron density distribution with changing 
R from R = to R = +00. Expression (jSJ together with 
a conventional assumption that a wave function and its 
derivatives are continuous functions gives us a solid basis 
for a consideration of such transformation. 

From the consideration given above follows that if $^, 
^B, and '^u are atomic wave functions then at i? = 
ca{0) = cb{0) = and cc/(0) — 1, while at R = -l-oo 
ca(+oo) = cb(+oo) = cn, where cn is a normalization 
coefficient and cu{+oo) = 0. Now, taking into account a 
continuous property of ^E'(-R), which means that all coef- 
ficients in (O and their derivatives are continuous func- 
tions on R, and suppose that a distribution of electron 
density is defined by a total wave function, we can derive 
two important statements: 

a) for any small Eq it is possible to find such Rq that 
when R < Rq the total molecular wave function can be 
presented as \I/(-R) = cu{R)'ifu + 0{eo)- This means that 
at sufficient small R the total molecular wave function 
and corresponding distribution of electron density must 
have only one maximum. As a consequence, a transfor- 
mation of the molecular wave function and the electron 
density with two maxima to those with one maximum 
passes through a flat wave function and corresponding 
flat electron distribution. 

b) the spherical component '^jj oi the total molecular 
wave function ^(i?) exists in molecule at i? 7^ 0. It is the 
wave function of a united atom and in molecule it repre- 
sents the electron density of a non-nuclear united atom 



or quasiatom. The spherical component of the molecu- 
lar wave function and appropriate electron distribution 
located at the middle of the bond can be observed in 
molecule at some R. 

Now, follow to Burrau [4|, we consider the Schrodinger 
equation in Born-Oppenheimer approximation for _ff^ 
in elliptic coordinates: £, = {ri+r2)/R and rj = 
[ri — r2) / R. In these coordinates, the wave func- 
tion for the E ground state can be presented as ^ = 
X(^, K)Y{ri, R), where X(^, R) and Y{t], R) satisfies the 
following equations: 
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and E and A is energy and separation constant to be de- 
termined. It needs to note here that the elliptic coordi- 
nates are undefined at i? = and therefore the equations 
given above are valid only for i? > 0. 

To clarify behavior of the total wave function at the 
middle point it is sufficient to investigate behavior of 
Y {f], R) and Y' (?/, R) a,t rj = and at small R. In this 
case equation Q transforms to an ordinary differential 
equation of the second order 
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The general solution of this equation is as a linear com- 
bination of two linear independent solutions: 

r(0, R) ^ ci (i?) exp (-VII77I) + C2 (i?) exp (\/^|77|) (4) 

that must be taken at 77 = 0. Despite of this limitation 
we can conclude that the algebraic structure of the to- 
tal wave function will be preserve at the vicinity of this 
point because of the wave function is continuous func- 
tion. This permits us to compare the expressions l@J and 
(0). From this comparison follows that the first solution 
of spherical symmetry corresponds to the wave function 
of a united atom which has a maximum on ry at 77 — 0, 
while the second one of cylindrical symmetry has no local 
minimum on 77 at ry = and is formed by the linear com- 
bination of wave functions of separated atoms. On this 
basis we can conclude that the both linear independent 
solutions are square integrable functions. Therefore, the 
general solution of the Schrodinger equation indeed is a 
linear combination of two linear independent solutions 
that represent two different quantum objects. 

It is interesting to note that Eq. Q) is similar to the 
one-dimension Schrodinger equation for an electron in a 
central potential 



dx^ 



{V -E)-^ ^0 



A formal solution of this equation is also a linear combi- 
nation of two linear independent solutions 



* = di exp -\/E-V\x\ + da exp y/E-V 



However, in this case the second solution is not square 
integrable function. Therefore, it is excluded from a con- 
sidered. 

Now, suppose that the coefficients in Q are continuous 
functions, we can conclude that at i? ^ Eq. Q can be 
smoothly transformed to the well-known one-dimensional 
atomic wave function only if C2{R) ^ (because of this 
component transforms to the not square integrable func- 
tion) and ci{R) becomes equal to a normalization con- 
stant of the united atom wave function. 

The first derivative of Q on 77 taken at 77 = is 

dYjO^R) 
drj 



= -VA[ci(i?)-C2(i?)] 



It is negative at small R because of C2[R) ^ at i? ^ 0. 
This means that Y (0, R) has a maximum. With increas- 
ing R and growing C2{R) the first derivative becomes 
equal to zero and after that it can be positive or neg- 
ative in dependence on the values of the coefficients. At 
large R when coefficient C2{R) becomes dominant only 
"minimum" can be observed at ry = 0. 

In the light of the analytical solution of the Schrodinger 
equation given above it needs to comment the results pre- 
sented in |y, |7| . Our analysis will base on the fact that 
the Schrodinger equation is a second order differential 
equation and its general solution is a linear combination 
of two linear independent solutions. An additional re- 
quirement for this equation is that the solutions must be 
square integrable functions. Otherwise they should be 
excluded from a consideration. 

An analysis of the publications [a, 13 shows that in 
these studies only one solution of cylindrical symmetry 
(that can be called classical) of the Schrodinger equation 
has been taken into consideration, while the existence 
of the second solution and its analytical properties was 
not discussed. This means that an implicit assumption 
that the second solution is a not square integrable func- 
tion has been used in these works. However, the present 
study gives strong evidences that the second solution of 
spherical symmetry (that can be called quantum or quasi- 
atomic) is a square integrable function and therefore it 
must be taken into consideration also. 

A proof that the results presented in [y, lZJ have been 
obtained only for the classical component of the total 
wave function but not for the total wave function fol- 
lows directly from the fact that the statement (v) of 
Theorem 1 [7| that the wave function has no local min- 
imum is a simple consequence of the analytical solution 
given above. 

To visualize the transformation of the total wave func- 
tion ^(i?), the electron densities of H2 at different R 



have been calculated in quantum mechanical variational 
calculations by using the Gaussian 98 program [3 . These 
calculations can be considered also as an independent 
verification of the analytical results presented above be- 
cause only the nuclear centered basis of 15s6p5(i4j3(;2/ilj 
Gaussian spherical functions was employed. It was 
formed from an optimized set of 15s functions supplied 
by polarization functions taken from pV7Z basis for H 
atom J3| • A quality of the variational wave function was 
controlled by comparing obtained variational total ener- 
gies with those ones calculated by the program [10| and 
presented in Table I. A uniform deviation of the vari- 
ational total energies from the numerical ones, which is 
equal to 1.2* 10~^ Eh, points out that the employed basis 
results in a high quality variational wave function. 

The calculated electron densities of H2 at a few R are 
presented on Fig. 1. At the beginning of a consideration 
we can note that at i? = the electron density is spher- 
ical. Then, at small i?, the spherical distribution con- 
tinuously transforms to an ellipsoidal distribution with 
one maximum. This distribution is presented on Fig. 1 
at i? = 0.008 a.u. Further increasing R results in appear- 
ing a flat surface on the electron density distribution at 
R = 0.010 a.u. given on Fig. 1 followed by appearing the 
two maxima at the positions of two nuclei and a transfor- 
mation of the ellipsoidal electron density distribution to 
a cylindrical one. Appropriate cylindrical electron den- 
sity at i? = 0.012 a.u. is presented at Fig. 1. Subsequent 
increasing R reveals the existence of the spherical com- 
ponent of electron density with additional maximum lo- 
cated at the middle of the bond. The three-dimension 
figure of the electron density with a quasiatom is pre- 
sented on Fig. 1 at i? = 0.019 a.u. 

To show that a quasiatomic solutions exist in any 
molecular system the electron density distribution has 
been investigated in all homonuclear diatomic molecules 
of the first-row elements. The spherical symmetry elec- 
tron density distributions (quasiatoms) arisen from the 
second linear independent solutions and located at the 
middle points have been observed in all molecules. At 
some internuclear separations quasiatomic solutions re- 
sult in additional maxima on a profile of electron densities 
along the molecular axis. An example of such electron 
density distribution in Li2 at 5.0 a.u. (i?e=5.051 a.u. 
[llj) is presented on Fig. 2. The quasiatom located at 
the middle point is explicitly recognizable on this figure. 
The electron density distribution was obtained in ah ini- 
tio calculations with cc-pVQZ Gaussian basis set [l3| by 
quadratic configuration interaction method as realized in 
Gaussian 98 program. Details of this study will be given 
elsewhere |l3| . 

The qualitative behavior of electron densities presented 
on these figures is fully consistent with that one, which 
follows from the analytical solution of equation (0J and a 
qualitative consideration of the asymptotic properties of 
the Hamiltonian (1) given above. Despite of using only 



nuclear centered basis, that corresponds to the cyHndri- 
cal solution, the spherical component of the molecular 
wave function or the other linear independent solution of 
the Schrodinger equation and the corresponding spherical 
electron distribution located at the middle of the bond 
was formed in accordance with theory of the second order 
differential equation to keep asymptotic and continuous 
properties of the total molecular wave function. 

Thus, the analytical solution of the Schrodinger equa- 
tion given in the present investigation shows that the 
ground state total wave function of H2 is a linear combi- 
nation of two linear independent wave functions of spher- 
ical and cylindrical symmetry. The former one repre- 
sents a quasiatom of the electron density of non-nuclear 
united atom and is responsible for appearing a spheri- 
cal distribution of the electron density and an additional 
maximum on a profile of the electron density along the 
internuclear axis at some R. 



TABLE I: The total energies (in Hartree) of the iifj molec- 
ular ion. 



R (a.u.) 


Etot (var.) 


Etot (numer.*) 


0.008 


123.00016807 


123.00016795 


0.010 


98.00026148 


98.00026136 


0.012 


81.33370830 


81.33370818 


0.019 


50.63250557 


50.63250545 



Calculated by the program pil | . 
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FIGURE CAPTIONS lar ion at different intermiclear distances. 

Figure 1: The total electron densities of H2 molecu- Figure 2: The total electron densities oi Li2 at 5.6 a.u. 
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Figure 1 . The total electron densities of //^"^ molecular ion at different intemuclear distances. 




Figure 2. The total electron densities of Z/2 at 5.0 a.u. internuclear separation. 



